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Abstract 

The Hubbard model is used to study an electronic system at half 
filling. Starting from a functional integral representation the spin- 
up Grassmann field is integrated out. It is shown that the resulting 
spinless fermion theory has an instantaneous cluster interaction, and 
that the spinless fermions are coupled to thermally fluctuating Ising 
spins. The coupling parameter of the spinless fermion interaction is 
a product of the Hubbard interaction and the hopping rate. As an 
example the strongly metallic as well as the strongly insulating regime 
are investigated in terms of the effective Ising statistics. 

PACS numbers: 71.10.-w, 71.10.Fd, 71.30.+h 
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I. INTRODUCTION 



The Hubbard model was originally constructed to describe a metal-insulator tran- 
sition for spin- dependent fermions in a simple way |T]-[§]. This transition reflects 
the competition between potential (static) energy and kinetic energy. The model 
is defined on a lattice, where the potential energy consists of a chemical potential 
and an on-site repulsion of fermions with opposite spin. The kinetic energy is given 
by a nearest neighbor hopping. It turned out from a number of calculations that 
this model has a rich structure because of the complicated interplay of charge and 
spin degrees of freedom. For instance, mean-field calculations for a magnetic order 
parameter indicate para-, ferro- and antiferromagnetic states for the half-filled sys- 
tem 0. Thus, the magnetic properties of the model became a central subject of 
investigations in solid state physics. 

The metal-insulator transition was discussed originally by Hubbard using self- 
consistent approximations Q, later in terms of a variational approach f|0], and in 
the limit of an infinite dimensional lattice |8||^,[|]. Very interesting investigations 
were obtained from computer simulations which indicate an insulating phase at 
half filling for sufficiently strong fermion interaction [10|,11|. However, the detailed 



mechanism and the properties of the transition are not entirely clear. 

To study the metal-insulator transition one can, in principle, start either from 
the metallic or from the insulating side. As the simplest approximations one could 
use non-interacting fermions on the metallic side or the local limit on the insulat- 
ing side, where the hopping rate is zero. Unfortunately, neither of these starting 
points is very useful in order to understand the interacting Hubbard model: Non- 
interacting fermions are unstable against an arbitrarily weak interaction , and the 
local limit is completely degenerate with respect to the spin. Therefore, an arbitrar- 
ily weak hopping rate would lift the degeneracy leading to a new state which might 



be magnetically ordered []T2j . The basic idea of the present work is to start from the 
extreme insulating as well as from the extreme metallic state at low temperatures 
and to construct a perturbation theory without analyzing its magnetic order. The 
latter is a restriction which simplifies the calculations significantly because the spin 
degree of freedom can be ignored. 

In this work a grand canonical ensemble is considered, where on average one 
fermion per site (half-filled system) is assumed. The non-interacting fermions as 
well as the static fermions (i.e. fermions without a hopping term) have a 2 M - 
degeneracy (M is the number of lattice sites) because each site can accomodate a 
fermion with spin-up or one with spin-down. Consequently, a perturbation theory 
around one of these limiting states is plagued by the degeneracies. For instance, a 
perturbation around the static state is a spontaneous hop of a fermion from any site 
to its nearest neighbor site. As a consequence, the fermion spontaneously creates 
a doubly occupied site and an empty site. The doubly occupied site may decay 
after some time again into two singly occupied sites. The resulting state is two-fold 
degenerate because of the possible two spin orientations. The unperturbed state can 
be an antiferromagnetic (Neel) state. A hopping process at a time ti can exchange 
two neighboring fermions which leads to two pairs of neighboring fermions with 
parallel spins. At time £2 the inverse hopping process can re-create the original 
antiferromagnetic state. Therefore, the two hopping processes are not independent. 
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Moreover, the intermediate state between time t± and t 2 has the same energy as the 
antiferromagnetic state. This implies a constant interaction in time. Consequently, 
the linked cluster theorem cannot be applied, since it works only for independent 
clusters or clusters which interact with a decaying interaction [|1^,|14]]. The central 
point of the present work is a concept which deals with this degeneracy. 

In order to control this exponential degeneracy it is natural to eliminate one 
spin orientation. This can be achieved formally by integrating out one of the spin 
orientation in the functional integral representation of the Hubbard model. The re- 
sult of this operation reveals an important structure of the effective spinless fermion 
model which is formally an expansion of the model in terms of the degeneracy: the 
expansion terms are not degenerate and the perturbation expansion can be applied 
independently to each of them. It turns out that the expansion is equivalent to the 
summation over the 2 M states of thermal Ising spins which are coupled to the spin- 
less fermions. After an approximation which is applicable for the strongly metallic 
and the strongly insulating regime of the Hubbard model, the fermionic degrees of 
freedom can be integrated out. Thus, the physics is described by the Ising spins: 
The strongly metallic regime is characterized by a ferromagnetic Ising structure 
in which the fermions can freely move at low temperatures. On the contrary, the 
strongly insulating regime is characterized by an antiferromagnetic Ising structure 
which creates a gap for the fermions, in formal analogy to a Peierls instability. 

The article is organized as follows: In Sect. II the Hubbard model is defined in a 
coherent state representation for a grand canonical ensemble of fermions. The static 
limit (no hopping) of the Hubbard is briefly discussed in Sect. III. Then in Sect. IV 
the integration over the spin-up component of the model is performed. The result- 
ing model of spinless fermions, which has a complicated but instantaneous cluster 
interaction, is analyzed in Sect. IV.A. In Sect. V the Ising spin representation of the 
spinless fermion model is introduced and discussed. Finally, in Sect. V.A the weak- 
coupling limit and in Sect. V.B the weak hopping limit are studied. Appendices A, 
B, and C give details of the calculations. 

II. THE HUBBARD MODEL 

The Hubbard model describes fermions with spin a =|, | on a ci-dimensional 
lattice A. It is defined by the Hamiltonian |2]|| 

H[ci{r), Cr7 {r)] = -t 4( r K(r') + 5I [f^Y, c l( r M r ) + Uc \( r ) c ^ r ) c \( r M r )}^ 

(r,r'),a r a 

where cj.(r), c CT (r) are fermion creation and annihilation operators, respectively. 
t > is the hopping rate, (r, r') means pairs of nearest neighbor sites on the lattice. 
fi is the chemical potential. 

Using this Hamiltonian a grand canonical ensemble of fermions at the inverse 
temperature j3 can be defined by the partition function, given in terms of a functional 
integral (coherent state representation) on a Grassmann algebra [jnj. For the latter 
the integration over a complex Grassmann field (^(r, t), ^o-( r ? t)) is given as a linear 
mapping from a Grassmann algebra to the complex numbers. At a space-time point 
(r, t) we have for integers k, I > 
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I [^ a (r,t)] k [^ a (r,t)] l d^ a (r,t)d^ a (r,t) = 5 M <5 M - 



The partition function of the grand canonical ensemble of fermions then reads 

Z = J exp(-S)T>[V,V] 

with the action 

S = tA E ^«(r, t) [tt ff (r, t) - * ff (r, * - A)] + zA £ ~#[*a(r, *), tt ff (r, t - A)] 

r,t ^ t h 

(1) 

and the product measure 

The discrete time is used with t = A, 2 A, /3. ^ ff (r, i) and ^(r, t) are independent 
Grassmann fields which satisfy antiperiodic boundary conditions in time ty a (r,(3 + 
A) = — ^ a (r, A) and ^ a (r, /? + A) = — ^ CT (r, A). For the subsequent calculations it 
is convenient to rename ty a (r,t) — > \fr a (r,t + A) because then the Grassmann field 
appears with the same time in the Hamiltonian of the action ([!]). 



III. THE LOCAL LIMIT 



Neglecting the hopping term in the Hamiltonian (i.e., for t — 0), the integration 
in the partition function factorizes in space, and the corresponding expression can 
be evalueated as 

Z = f nexph^Mp^J = I] / IIexp[-^M)P[^J = Z? 1 , (2) 

r,t r J t 

where 

Sa = E {^( r > t)**(r, t+A)-fl^ a (r, t)V a (r, t)+iAU^(r, t)tt T (r, t^^r, t)^(r, t)} 

t,r 

with h — 1 and p, — 1 — zA/x. Z\ is the partition function of the Hubbard model 
with one lattice site (static, local or atomic limit): 

Zi= f Y[e^[-S d {r,t)]dti{r,t) = 1 + 2^ A + [ft 2 - iAUf /A . 
•* t 

Using the new parameters // = iA/j,, U' = iAU, and /?' = /3/A, we can define the 
following weights, depending on the number of particles per site 

w = l/Zi 

Wl = 2(1 -fi'f/Zt 

w 2 = [(l- f i') 2 -U'f/Z 1 . 
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In the temperature formalism, where the time t is replaced by the imaginary time 
through a Wick rotation, the weights w , W\ and w 2 are statistical weights. Then 
the average number of particles per site is 

n = Wi + 2w 2 - 

At zero temperature {(5' — > oo) this gives 

f if < /i' < 1 + VTTu 7 
n = < 1 if fx' < and //(// - 1) < U' ■ 
{ 2 otherwise 

The groundstate of the expansion around the local limit is degenerated because 
each singly occupied site can accomodate a fermion with either spin j or with spin 
j. This degeneracy must be handled with care. In particular, to obtain a unique 
hopping expansion one has to separate the degenerate contributions. Therefore, one 
cannot directly work with the hopping term as a perturbation but have to set up 
a perturbation theory which works with non-degenerate groundstates. This means 
that one has to divide the degenerate groundstate such that the perturbations remain 
in their corresponding non-degenerate groundstates. 



IV. INTEGRATION OVER THE SPIN-UP FIELD T 

The action S can be divided into three pieces as 

S = 5*| + Si + Sj 

with 

S ° = E { E [*-(r, t)**(r, t + A) — ji* a (r, t)* ff (r, t)\ - r £ * ff (r, t)V a (r', t)} 

for a — t, J, with r = iAt. The interaction between the two spin orientations is 
given by 

Sj = f/'5:^ T (r,t)M/ T (r,t)^(r,t)^(r,t). 

r,t 

Now it is possible to integrate out the spin- up field y&f, since the field ^| appears 
in S only as a quadratic form. The integration over this Grassmann field gives a 
determinant 

f e~ s i~ Sl Y[dVi(r,t)dVi(r,t) = det[- d t + p, + i - U'Vl ^J, (3) 

J r,t 

where dt is the time-shift operator 



t = p 

The last equation is due to the antiperiodic boundary condition of the Grassmann 
field. This definition gives 
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(dt)- 1 = cf , det(dt) = I. 

We assume that the number of time slices f3' is even such that det(—d t ) = det(d t ) = 1. 
The matrix t T y = r if r, r' are nearest neighbors and zero otherwise. Expressions 
in the determinant which do not have a specified matrix structure are implicitly 
multiplied by the corresponding unit matrix. For instance, p, is multiplied by the 
space-time unit matrix whereas t is multiplied by the time-like unit matrix. 

In the following subsets of the space-time lattice A = A <E> {A, 2A, f3} will be 
considered. For a subset A k C A we define the determinant of the the projected 
matrix P k AP k as 

det Ak A = det Ak (P k AP k ), 
where P k is the projector onto A k . 



A. Effective Cluster Action of Spinless Fermions 

The partition function is now a functional integral of the spin-down Grassmann 
field 

Z = J e' s ^det[ + £ + U'V^dV^A. 

Formally, the determinant could be expressed as part of the action by using the 
identity detA = exp[Tr(log A)]. However, this would be too naive because the term 



Tr 



log(- dt + fi + i-U'V^L) 



has a complicated interaction of the Grassmann field in space and time. Moreover, 
at least for t — the interaction has a long-range part in time which reflects the 
degeneracy of the unperturbed system. Fortunately, there is a way to avoid these 
difficulties: As shown in Appendix A, the determinant can be expanded in terms of 
the partitions A k C A of the space-time lattice A as 

detl-dt+fi+i-U'V^d = J2 det A J-(/Z+f)^]exp[rr Afc log(l-(f+/Z)- 1 i7 , * i *;)' 

A fe CA 

The partitions include the empty set which gives det^A = 1. This expansion is the 
most important step for the treatment of the Hubbard model in this work. The first 
consequence is that the partition function Z is now given by an expansion in terms 
of A k as Z = E Afc Z Ak with 

Z Ak = det Ak [-(jl + t)dj] J exp(-S A jZW ( 4 ) 

and the action 

S Ak = * A • (d t - fl - f> A - Tr Ak log (l - (t + ji)- 1 !/'*^). 

The second term of S Ak can be expanded in powers of the Grassmann field. This 
yields an instantaneous cluster interaction 
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5i = Ej ^ E (^^)^^l(n,t)^(ri,t)...(U/i)-i iir ^(r,t)^(r,t). 
J>1 (r,t)eA fe n,...,n-i 



Due to the identity 



(i+p)- 1 = -- -{i + Ji)- 1 



the action 5j can also be written with u = U'(i + /i)^ 1 as 

(r,i)GA fc 



E-T E E (-(t + ^)- 1 ) rri ^(r 1 ,t)vI/ i (r 1) t) 

l>2 1 (r,t)&A k ri,...,n_! A* 



...(-^(U/i)- 1 )^*;^*)^^*). 

This shows that the interaction terms with I > 2 do not distinguish between weak 
interaction (If ~ 0) and weak hopping (t ~ 0) because the expansion parameter is 
U't. 

The non-interacting limit U' = as well as the local limit t = can be checked 
immediately because of a vanishing interaction. The former gives the product of 
two determinants (from both spin orientations) 



Z = det{ -d t + p + iy 



and the latter 



Z = Z[ 



M 



in agreement with Eq. @). 

V. ISING SPIN STATISTICS 

It is interesting to notice that the determinant det\ k [— (ft + t)df] in Z\ k deter- 
mines which partition A& contributes with non-zero weight. Since dt is diagonal 
in space but off-diagonal in time, the determinant det\ k (df) is non-zero only if a 
site r is in A& at all times. Therefore, the contributing partitions are of the form 
Afe = \k x { A, 2A, ...,/?}, where A^ is a partition of the space lattice A. The deter- 
minant then reads 

det Ak [-{fl + t)dJ] = [det Xk {fl + t)f. 

Now it is convenient to introduce the projector I\ k which projects on A& 

it \ fl ifreA fc 
1 Ajr 10 ifr^A fc ' 
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Consequently, the determinant can be expressed as 

[det Xk (jl + t)f = [det x (l + I Xk [ft + t - l]I Xk )f- (5) 
The interaction can also be written in terms of the projection operator: 

1>1 t r,n,...,rj-i 

The partitions A/, enter Z only through the projection operator I Xk . Eqs. (|J) and 
(D imply 

2 = E ^ = E [**a(1 + h k [fi + t - ^jf I expC-^jP^] 

A fc A fc 

with the action 

5 Afc = ^ • (ft - t - ft + uI Xk )* l + 0{{U ! i/p?Y). 

The projection operator J Afc can be expressed as (I({S r }))„j = (1 + S r )5 rr >/2 with 
the Ising spin 



S <■ 



1 r E Afc 
-1 r £ Afc 



The partitions Afc enter Z only through the projection operator I Xk . Therefore, 
the sum over the randomly chosen partitions Afc of the space lattice is equivalent 
to the sum over randomly chosen Ising spins. This is a thermal Ising model with 
complicated interaction spin interaction in a magnetic field. The latter is implied 
by the missing invariance under a global transformation S r — > —S r . 

Neglecting the interaction Si, which is weak for the metallic as well as for the 
insulating regime, one gets 

Z^Z I= [det x {l + + det{-d t + ji-ul + t). (6) 

{S r =±l} 

The partition function Zj can be used as a starting point for an approximative 
treatment of the Hubbard model. In the rest of this work the favoured spin struc- 
tures will be analyzed, i.e., the spin configurations with maximal Boltzmann weight. 
Details of the calculations are given in App. B and C. 



A. Weak Interaction: U' ~ 

A metallic behavior is expected in this regime. Since 

[det Xk (fi + t)f 

is independent of U', only the second determinant of Zj must be expanded around 
U' = 0. This yields in leading order 
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det(-d t + p-ul + i)^ exp[ / 0(|w| - 1) \og(n)d d k] exp^'/i^l + S r )]. 

J r 

with the effective magnetic field 

u r Q(\ K \ — I) d 

h = - / — -d d k with k = p + 2t ^ cos kj. 

2 J K 

Thus, the weakly interacting regime prefers a ferromagnetic state in terms of the 
Ising spin S r — 1. Knowing the Ising spin configuration with the highest weight, 
one can return to the partition function of Eq. ([]) and obtains 

Z l ps [det(# + £)] /3 'det(-9 i + p - u + i). 

The argument of the determinant can be considered as an inverse Green's function 

G = (-d t + p, - u + t)- 1 . 

Thus, the effect of the ferromagnetic Ising spin structure is a constant shift of the 
chemical potential. Of course, the fluctuations of the Ising spins around the structure 
with maximal weight have to be taken into account as soon as j3' < oo. 



B. Hopping Expansion: t m 

The expansion of the first determinant of Zi in powers of r is 



^/3'|A fc | exp [ 



w (t; "" )2] 

r,r€\ k 



(7) 



and the expansion of the second determinant is 



det(—d t + fi — ul + t) p» n 



-/3'(M-|A|) 



exp 



/3V 2 



(r,r) 



Here the potential terms V r , V T > are either jt/or/i- u ~ fi — U' '/ ' \i with the function 

W(V r ,V r ,) -- 



for V r = V r i = fi — u 

-1 1 (up) for V r ^ V r > 

1/p 2 for V r = V r r=p 



For r = all terms in the sum J2\ Z\ have the same weight. This reflects the fact 
that there is a degenerate perturbation theory for the hopping expansion in the 
half-filled system due to the spin degree of freedom. For r/0 the second factor of 
(H) can be approximated for (3 ps oo (j3 is now real) as 



exp 



exp 



(r,r) 



4p* 

(r,r) 



(9) 



According to (0) and @ the maximal contribution to Zj comes from V r ^ V r > on 
nearest-neighbor sites. This corresponds with an antiferromagnet state of the Ising 
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spin S r = (— l) ri + - + r d. The partition function in Eq. @ reads with the staggered 
Ising spin structure 

Z, « 2^' M ' 2 det{-d t + V AFM + £), 
where Vafm is the staggered antiferromagnetic potential 

v AFM = fi-^[i + (-ir + - + ^}. 

As an effective Green's function of the fermions one can study 

G=(-d t + V AFM + i)-\ 

The staggered potential creates a gap. This can be seen, for instance, in the special 
case of a two-dimensional lattice. Then the eigenvalues of G~ l are 



u ~ 2 ± VT + |r/l1212 

with h\2 = l + e lkx +e lhy -\- e lk *+ lk y [—^ < fc. < yj-). This result describes the physics 
of an insulator with a gap 

U/2, 

in contrast to the homogeneous potential of the ferromagnetic Ising structure for 
the weak interacting case. 

There can be a phase transition from weak coupling (ferromagnetic Ising system) 
to weak hopping (antiferromagnetic Ising system) even at finite (3 (i.e., non-zero tem- 
perature). It is not clear, however, whether this is a first or second order transition. 
Numerical (Monte Carlo) simulations may reveal more details of the transition. 

VI. CONCLUSIONS 

Considering quantities of the Hubbard model which depend only on one spin 
direction, e.g., spin-down, it is possible to integrate out the other spin direction 
in the functional integral. This idea was carried out for the partition function Z. 
The latter could have been extended to a generating function for spin-down Green's 
function without modification of the procedure. The result of this integration is 
a functional integral which has a representation in terms of thermal Ising spins. 
Moreover, the fermionic interaction of the effective model has a coupling parameter 
U'i which becomes small in the limit of weak Hubbard interaction U' as well as in the 
limit of low mobility of the particle (small hopping rate t). Therefore, the effective 
fermionic interaction can be neglected if one is only interested in the strongly metallic 
and the strongly insulating regime at half filling. The advantage of this approach 
is that the degeneracy of the special cases U' — and t = 0, which are difficult 
to handle, are controlled by the coupling of U' and t to the Ising spins: a small 
Hubbard interaction U' creates an effective magnetic field h oc U' which couples 
linearly to the Ising spin. On the other hand, a small hopping term favours an 
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antiferromagnetic (staggered) Ising spin configuration because of the effective Ising 
spin interaction 



t 2 

\r,r) 

The Boltzmann weights of the effective Ising spin system were used to extract the 
most relevant contributions to the sum over the spin configurations. 

The introduction of Ising spins by a direct decoupling of the interaction term in 
Eq. ([!]) (Hubbard-Stratonovich transformation) is a well-known approach. It leads 
to a dynamic Ising spin . The approach of the present work avoids the dynamics 



of the Ising spin. However, the price for this simplification is a cluster interaction 
of spinless fermions. 

Although the approximations used in this work do not provide insight into the 
transition from the metallic to the insulating phase, a more accurate treatment of 
the partition function Zf of ([]) by including fluctuations of the Ising spins, e.g., 
using a Monte Carlo simulation, may allow to access the transition at not too low 
temperatures. 
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APPENDIX A 



The space-time determinant on the r.h.s. of Eq. @ can also be written as 

det[ -dt+fL + t- U'V^d = det(l + A) = £(-1)* J] |X(r,t),(r,t) + A-,*r(r,t)j 

(r,i)GA 

with the matrix A = —(p, + t — U'^ ^ t)df . The product over the lattice sites gives 
a sum over all subsets A& C A of the space A and their complements A' fe = A \ A& 

Yl XX _i r[ II A r,t;n(r,t) [ Y[ K(r,i),{r,t) ■ 
AfcCA tt (r,t)eA k (r,t)£A' k 

The Kronecker delta S^^^t) on A^ implies ir(r, t) G A k for (r,t) G A&. Therefore, 
only that part of the matrix A contributes which is projected onto A&. This implies 
an expansion of the determinant in terms of all partitions of the space-time lattice 
A as 

det[ -d t + fi + i- U'*^] = det Ak {P k AP k ) = ^ det Ak A (10) 

A fc CA A fc CA 

with det$A = 1 for an empty set A&. The projected determinant is 
det Ak ( - (fl + 1 - U'V^Jdj) = det Ak [-{fl + t)dJ]det Ak (l - (fi + t)~ l U'^^ l 
where the second second determinant reads 

exp [Tr Ak log (l - (t + fLy'U^^^ . 



I) i 
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APPENDIX B 

Expanding the second determinant of Zi up to first order in 5V yields 
det(-d t + jl + 5V + i)^ det(-d t + fl + t) exp (Tr[6V(-d t + Ji + i)' 1 
The trace term can be written as 

r 

where 

u f G(\k\ - 1) d d k . , n * 

h — — — -j with k = a + 2r > cos kj. 

2 J k (2n) d P ^ 3 

APPENDIX C 

The first determinant of Zi can be expanded up to second order in t as 
[det Xk (fl + i)f «^lexp[-^ ]T (^') 2 1- 

The second determinant of can be expanded in powers of t up to second order 

det(-d t + ft + 5V + i) 

« det(— Sfe + /Z + <jy) exp ( - ^Tr[i(-d t + Ji + 5V)- l i{-d t + p + 5V)- 1 ]) . 
The trace term reads 

Tr[i(-d t + fi)-H(-d t + fi)- 1 ] = (3>Y,(irr>) 2 r t^— 

where the integral gives 

f 2 * 1 1 rf w 



1 duj 



r e w V r > - e w 2n 



o K - e iu} V r > - e iu) 2ir 



' for V r = V r , =\x- U'/jl 

— l/(up) ioiV r ^V r i 
k l//i 2 for V r = \4' = At 



Due to r 0, the approximation w t^'/A 4 can be used. Then for f3' m oo the 
determinant is 

det(-d t + ft + 5V)= + + SV r f) » ^'( M -I^D, 

r 

since |/2| > 1 and |/i — {7'//i| < 1- 
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